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The strong coupling constant o s (/^o), taken at a fixed reference scale fio, is the 
single free parameter of QCD and should be known to the highest available preci- 
sion. The value of a s should also be determined with good accuracy over as large 
a range of scales as possible, in order to reveal potential anomalous running in the 
strength of the strong interaction. Lattice QCD is now able to calculate a s with 
accuracy comparable to or better than experiment. We review the status of such 
lattice calculations in quenched and full QCD. 



1 QCD coupling and /^-function 



The scale dependence of the QCD coupling a s 
/3-function: 
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where Pq — 11 — 27V//3, pi = 51 — 19./V//3 and Nf is the number of flavors of 
quarks with mass less than the energy scale [i . The coefficients /3 n for n > 2 
are renormalization-scheme-dependent . 

By solving the differential equation for a s , a constant of integration is 
introduced. This constant is the one fundamental constant of QCD. Usu- 
ally one chooses for this constant the value of a s at a fixed reference scale 
fio. Equivalently one can consider the dimensional parameter A, defined for 
example as 
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Note that, by using this relation at the leading order, one gets that A(JV/) is 
known with an accuracy of about 5% if a s (/z) is given with an accuracy of 
about 1%. 
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When two different rcnormalization schemes are considered one can write 



and check that 
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Using the A parameter we can also obtain a parametrization of the fi depen- 
dence of the QCD coupling" 
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Clearly the coupling decreases as the scale \x increases, corresponding to the 
property of asymptotic freedom. However, due to the logarithmic dependence 
one needs to consider a large range of scales \i in order to see a variation of 
a s {n) by a factor two or three (see Fig. 9.2 in Ref. 2 ). 



1.1 Experimental results for a s (/z) 

The value for a quantity A from an experiment can be used to extract a 
value of a s {n) by considering the following perturbative expansion (in some 
renormalization scheme) 

A = c + cia s (ji) + c 2 a 2 s (p) + . . . . (7) 

If the scale involved is low one should also try to take into account non- 
pertubative contributions to A . Clearly the physical result A does not de- 
pend on the chosen renormalization scheme. However, the truncated perturba- 
tive scries does exhibit renormalization-scheme dependence. In particular, the 
finite coefficients Cj (i > 2) depend on the renormalization scale /j, and, im- 
plicitly, on the chosen renormalization scheme. Therefore, for a given scheme, 
one should find the best choice 6 for the scale /j, in order to reduce the scale 
ambiguity in the expansion (7). 

In order to compare the values of a s (fj,) from various experiments 3,4 , these 
values must be evolved — using the renormalization group — to a common 

"Note that this parametrization is not unique but depends on the choice of the constant of 

integration (see for example problem 13.11 in Ref. x ). 

b See Ref. 3 for various methods proposed for choosing the scale ^i. 
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scale (usually the mass of the Z° boson) and related to a common scheme 
(usually the MS scheme) . The present world-wide average 3 in the MS scheme 
for a s (M z ) (Nf = 5) is 0.1172 ± 0.002. 

2 Numerical evaluation of a s (fi = q) 

One can perform the experiment also on a computer, using Lattice QCD. 5 For 
example, one can evaluate numerically the plaquette Wi t i = < Tr?7 M „/3 > 
— which corresponds to the elementary gluonic action on the lattice — and 
compare the result with a perturbative expansion. Also in this case one should 
be careful with the choice of the scheme used for the perturbative expansion 
and with the choice of the scale /x. For example, it is well known 6 that the 
expansion of the plaquette Wi ; i in terms of the lattice bare coupling a = 
5q/47t yields a very poor result. In fact, for g 2 = 1 one- loop perturbation 
theory with Nf — predicts — logWi.,1 « 0.4228 while the corresponding 
Monte Carlo result is about 0.5214. 

A better expansion 6 is obtained if one considers an effective coupling 
defined in terms of a physical quantity. For example, if V(q) is the inter- 
quark potential one can define 

V(q) = (8) 

and work in the ay scheme. In this scheme one gets, at the one-loop level 
and with Nf = 0, 

-]QgWi,i = c iav {q 2 ) [1 + c 2 a 2 v {q 2 ) + ...], (9) 

where c 2 = log (^r^j + d and d « 1.329 . In order to fix the value of the 
scale q the idea introduced in Ref. 6 is to impose the relation 

a V (f) Cl - a v (f) J d\f(q) = J d\a v {q 2 )f{q) , (10) 

namely q is the average momentum entering in the one-loop calculation. By 
expanding both ay(g 2 ) and ay(q 2 ) in terms of ay(/i 2 ) one gets 

log(g 2 ) J d 4 qf(q) = J d 4 q log (q 2 ) f(q) . (11) 

In the case of the plaquette this prescription gives q w 3.4018/a where a 
is the lattice spacing. Clearly, this procedure can be repeated for different 
Wilson loops W r> t, i.e. from the numerical values for W r ^ one can obtain 
av{q 2 t ) using a perturbative expansion like eq. (9). As can be seen from Fig. 
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8 in Ref. 6 , the coupling ay(g 2 ) obtained from such calculations runs with q 
in agreement with a two-loop formula (6). 

In order to compare these results with experimental data one needs to 
fix the lattice spacing a, i.e. the scale q should be given in physical units. 
This can be done by comparing lattice data with an experimental input, as 
has been done 7 using the T meson spectrum, known experimentally with very 
good accuracy and for which accurate simulations can be done. Using results 
for Nf = and Nf — 2 dynamical fermions, which allows extrapolation to 
physical case Nf = 3, the authors of Ref. 7 obtained a final result — in the MS 
scheme and evolved perturbatively to the Z° mass scale — of 0.1174 ±0.0024. 
Similar results have been obtained by other groups, using different lattice 
setups, and have been included by the Particle Data Group 3 in the evaluation 
of the world-wide average of a s {Mz) <<N,= ^ ■ 

2.1 Ingredients for a numerical evaluation of the QCD running coupling 



From the previous analysis it is clear that in order to evaluate the QCD 
running coupling constant using lattice numerical simulations one needs: c 

• a good definition for a s allowing high accuracy, control over finite-size 
effects and discretization effects, 

• an accurate determination of the energy scale (i.e. the lattice spacing a ), 

• a perturbative relation with ct-jjg, 

• the possibility of extending the procedure to dynamical fermions (un- 
quenched simulations). 

3 Running coupling in a finite box 

The Alpha collaboration 9 studies the QCD running coupling using finite boxes 
(with physical size L = Na), a definition for a s that explicitly depends on 
the scale L, and a non-perturbative finite-size scaling technique. 10 To this end 
one could consider, for example, the force between static quarks separated by 
a distance r = L/2, given by 



constant 



a q q(L) 
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r=L/2 



(12) 



c See Ref. 8 for a previous review of the numerical determination of a. 
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or a correlation function 11 for a separation r = L/2 . Most of the simulations 
by the Alpha collaboration have been carried out using as definition for a s 
the response of the system to a change of the boundary conditions in the 
Schrodinger functional approach (for details see Rcf. 12 ). In all cases the 
physical size L of the lattice is the only physical scale on which these couplings 
depend, i.e. the box size may be regarded as a reference scale, similar to the 
renormalization scale ji . Thus, one has a coupling running with /i = 1/L. 

The key point in considering these couplings is the use of a non- 
perturbative renormalization group technique. More specifically, one defines 
the step-scaling function a(u) through the relation a s (2L) — o~(a s (L)). 
Clearly this function provides the value of the coupling in a box of size 2L 
given a value for the same coupling in a box of size L. Thus, <r(u) can be 
regarded as an "integrated" /3-function. The step-scaling function can be eval- 
uated non-perturbatively using numerical simulations by considering several 
values of L, as shown in Fig. 4 in Ref. 13 . Notice that for each data point in 
the plot one should carefully check for possible systematic effects, especially 
discretization errors. In fact, each point is the result of an extrapolation^ to 
zero lattice spacing of data obtained from several simulations, done with the 
same physical size L but different lattice spacing a and number of sites per 
direction N. 

Once the step-scaling function is known, i.e. after fitting the numerical 
data, one can evaluate the running coupling in a very large box Lq = L max 
and then find the value of the coupling for boxes of sizes Li = L /2 l , using 
a recursive procedure and a numerical inversion of the step-scaling function. 
[In order to keep errors under control one should avoid using the fit for the 
step-scaling function outside the range of the numerical data, i.e. one should 
avoid extrapolating the step-scaling function a(u) to very small values of u] 
In this way one obtains the values a s (ft) for the running coupling at the 
scales fa = 1 1 Li = 2 l /L max . 

Since all the scales fa are given in term of L max , the last step consists in 
setting the scale for the largest lattice, i.e. finding the value of L max in phys- 
ical units by comparing lattice data with an experimental input. Thus, the 
use of a non-perturbative finite-size scaling technique avoids the requirement 
of a non-perturbative evaluation of the scale Li for very small values of Lf. 
this is in general complicated by finite-size effects and by the fact that small 
values of Li correspond to large value of (3 and that simulations at large /? 
and small physical volumes are essentially "perturbative" . 

The final value 13 obtained by the Alpha collaboration in pure QCD is 



d See Ref. 14 for possible problems with this extrapolation. 
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Ai^ = 251 ±21 MeV. Preliminary data with dynamical fermions (Nf = 2) 
have been presented last year 15 . 

4 Numerical determination of A 

Several groups study the running of the QCD coupling (in some scheme) 
over a range of energies and then try to fit the data using a two- or three- 
loop expression for a(q) as in eq. (6), considering A as a fitting parameter. 
Alternatively, one can evaluate A as a function of a(q) and look for a plateau 
at high scale [note that in eq. (2) the left-hand side is constant, for a given 
value of Nf, even though the right-hand side depends explicitly on the scale 

4 

4-1 Inter-quark potential 

One can evaluate the QCD coupling using the static potential V(r), which 
can be evaluated for several values of r with good accuracy using Wilson loop 
data and the relation 5 

W{r,t) = C(r)e [ -* y(r)1 + ... . (13) 

This formula should in general be considered in the limit of large t, but using 6 
smearing and noise-reduction techniques one can increase the overlap of the 
evaluated quantity with the ground state, i.e. increase the value of C(r), and 
obtain a good signal already for small values of t. Then, from V(r) one can 
obtain the force F(r) [and the running coupling a q q(l/r), see eq. (12)] using a 
finite-difference approximation for the derivative dV(r)/dr. In the quenched 

case 17 one obtains A^- -0 ^ = 246±7±3 MeV, in excellent agreement with the 

MS ° 

result from the Alpha collaboration. Unquenched simulations (also interesting 
to study string breaking 18 ) are under way 19 . 

4-2 Three- gluon vertex 

In this case one evaluates the three-gluon vertex 20 in the MOM scheme (in 
Landau gauge). For this quantity perturbative results are known to three 
loops 21 (and the corresponding beta function to four loops). In order to obtain 
good data a careful analysis of breaking of rotational symmetry is required 22 . 
From quenched simulations 23 one obtains the value A^—L~ n ^ = 295 ± 5 ± 15 

MS 

MeV, somewhat larger than the previous results. Preliminary unquenched 

e For a review of this method see for example Ref. 16 . 
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data 24 have been presented last year and the authors quote A-jJL = 264±27 
MeV. A similar study has started recently using the quark-gluon vertex 25 
yielding A^ =0) = 300 ±J|g ± 55 ± 30 MeV, in agreement with the result 
obtained using the three-gluon vertex. 



4-3 Gluon propagator 

In this case 26 one evaluates the gluon propagator D(q) (in Landau gauge) and 
defines Z 3 (q) = q 2 D(q). Then, the gluon data can be fitted by considering 
the two coupled differential equations 
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where T(a) is the anomalous dimension of the gluon renormalization constant 
in the MOM scheme. At the lowest order these equations imply the well- 
known result D(q 2 ) ~ l/q 2 log c (g 2 /A 2 ) with c = ^ = 13/22. 

A solution for these two equations depends on the values Z 3 (n) and 
a(n) given at some scale (1, and a direct fit of the data provides a value 
for these two constants. The value a s (n) obtained from the fit can then be 
evolved, using the renormalization group, to the Z° scale and related to the 
MS scheme, using the perturbative relation between the two schemes. The 

authors 26 quote A^ =0) = 319 ± 14 ±20 MeV as nnal result in thc quenched 
case, in agreement with the results quoted in Sec 4.2. 



4-4 Search for power corrections 

Recently there have been various attempts of finding power corrections to the 
perturbative running of a s , following theoretical predictions 27 . In particular 
data from the inter-quark potential 28 , the three-gluon vertex 29 ' 30 and the 
gluon propagator 30 have been reanalyzed. In all cases the authors found that 
a nonzero contribution proportional to 1 /p 2 improves the fit of the data. In 
particular, the result obtained in this way for the A parameter using the three- 
gluon vertex 29 , i.e. A^^. -0 ^ = 237 ± 3 °_ 10 MeV, is in much better agreement 
with the results quoted in Sections 3 and 4.1 than the original value (see 
Section 4.2). 
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5 Running couplings at the IFSC-USP 



We now describe the simulations being performed since July 2001 at the In- 
stitute of Physics of the University of Sao Paulo, Sao Carlos (IFSC-USP). 
In connection with a grant from FAPESP ("Projeto Jovem Pesquisador" , to- 
gether with Tereza Mendes), we have set up a dedicated PC cluster at the 
IFSC-USP. The system has 16 nodes with 866 MHz Pentium III CPU and 256 
MB RAM memory (working at 133 MHz), a server with 866 MHz Pentium 
III CPU and 512 MB RAM memory (working at 133 MHz) and is operating 
with Debian Linux. The machines are connected with a 100 Mbps full-duplex 
network and in total there are 130 GB in /tmp directories and 24.8 GB for 
the /home directory. 

We are carrying out simulations considering two definitions of the running 
coupling constant. In the Coulomb gauged we consider 31 

9 2 C (P) = §PMP) = ^.9o 2 P 2 ^44( P ) , (16) 

where V(p) is the color-Coulomb potential, D^p) is the time-time com- 
ponent (at equal time) of the gluon propagator and the momenta are three- 
dimensional. 9 Clearly, if V(p) is governed by a string tension at large dis- 
tances, i.e. goes like 1/p 4 at small momenta, then we should find g^{p) ~ 
1/p 2 in the infrared limit. In Landau gauge 71 we consider 32 

gl( P ) ee |L [p*D(p)] [p*G( P )] 2 , (17) 

where D{p) is the transverse gluon propagator and G(p) is the ghost propaga- 
tor. This expression enters directly in the quark Schwinger-Dyson equation 33 
and can be interpreted as an effective interaction strength between quarks. 2 
In this case one obtains g\{p) ~ p~ 2 if, for example, D(p) ~ const and 
G(p) ~ p in the infrared limit. On the contrary, if the gluon propaga- 
tor goes to in the infrared limit and the ghost propagator blows up not 
faster than p~ 4 then g\{p) has an infrared fixed point. The last behav- 
ior has also been obtained 32 by solving (approximatively) a coupled set of 
Schwinger-Dyson equations. 

•^Work in collaboration with D. Zwanzigcr, New York University (NYU). Some of our 
simulations are done on a PC cluster at NYU. 

9 Remarkably, this definition of the running coupling constant is renormalization-group 
invariant . 

h Work in collaboration with K. Langfcld and J. C. R. Bloch, Tubingen University. 
'This is also a renormalization-group-invariant quantity since (in Landau gauge) Z\ = 1. 
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In our project we are also interested in a numerical verification of Gribov's 
confinement scenarios — written in terms of the infra-red bahavior for the 
propagators considered above — for these two gauges. Theoretical studies 34,35 
predict that in Landau gauge, there should be a strong suppression of the 
(unrenormalized) transverse gluon propagator D(p) in the infrared limit and 
an enhancement of the ghost propagator G(p) in the same limit. These results 
clearly indicate the absence of a massless gluon from the physical spectrum 
and provide an indication of a long-range effect in the theory that may result 
in color confinement. The confinement scenario is particularly simple in the 
minimal Coulomb gauge where the ghost propagator determines directly the 
Coulomb interaction 34,31 . In fact, in this case, confinement of color, i.e. the 
enhancement at long range of the color-Coulomb potential V(R), is due to 
the enhancement of ghost propagator G(p) at small momenta. At the same 
time, the disappearance of gluons from the physical spectrum is manifested 
by the suppression at p = of the propagator Dij(p,p4) of 3-dimensionally 
transverse would-be physical gluons. 

Preliminary results in the pure SU(2) case for the two running coupling 
constants and Gribov's confinement scenarios have been presented in Refs. 36 . 
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